Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1069 1566 (N36E05 L= 76)   785    239        3 (1020/1451-1099/1643) -53.1 E10-W00 N28-N41
      2  714 1539 (N29E30 L= 52)   808    239        1 ( 656/1463- 793/1615) -31.0 E36-E21 N24-N33
      3  411 1430 (N17E49 L= 32)  2558    518        6 ( 367/1257- 460/1575) -45.0 E59-E40 N07-N26
      4 1710 1254 (N23W46 L=127)  1149    255        2 (1639/1134-1778/1364) -47.3 W39-W52 N16-N30
      5 1368 1241 (N18W19 L=101)   255     74        0 (1354/1205-1386/1286) -78.7 W17-W20 N15-N20
      6 1094 1228 (N14W02 L= 83)   202     68        1 (1053/1216-1134/1238)  14.0 W00-W04 N13-N14
      7  682 1159 (N04E23 L= 59)   216     74        0 ( 644/1147- 723/1170)  21.8 E25-E20 N03-N05
      8 1639 1109 (N13W40 L=121)   424    109        2 (1623/1059-1653/1170)  90.0 W38-W41 N10-N17
      9 1728 1103 (N14W47 L=129)   447     79        1 (1706/1067-1748/1140)  78.7 W45-W49 N11-N16
     10  195 1007 (S11E61 L= 20)  1215    504        2 ( 101/ 980- 306/1030)   6.1 E77-E49 S13-S08
     11 1702  655 (S14W54 L=135)   517    117        1 (1658/ 622-1734/ 671)   0.0 W49-W57 S15-S12
     12 1143  425 (S37W18 L=100)   915    173        4 (1050/ 401-1224/ 449)  17.1 W11-W24 S38-S34
     13  462  407 (S50E42 L= 40)   480    134        0 ( 436/ 380- 498/ 428)  26.6 E46-E36 S53-S47
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y