Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1752 1478 (N43W45 L=143)  1230    188        2 (1665/1442-1829/1509)  13.0 W33-W59 N40-N47
      2  372 1239 (S09E47 L= 52)   878    215        4 ( 300/1181- 446/1307)  48.0 E51-E42 S13-S03
      3  642 1272 (S02E28 L= 71)   219     73        1 ( 619/1244- 669/1305)  90.0 E29-E26 S03-N00
      4  264 1205 (S12E56 L= 43)  1150    161        5 ( 211/1133- 320/1264)  41.6 E60-E51 S16-S08
      5 1813 1025 (N15W52 L=151)  1264    257        2 (1785/ 905-1839/1140) -59.7 W47-W55 N08-N20
      6 1885 1013 (N17W61 L=160)  1208    220        3 (1862/ 927-1932/1123) -28.1 W56-W71 N12-N23
      7  266  908 (S30E52 L= 47)   721    237        2 ( 210/ 820- 325/1008)  45.0 E60-E46 S36-S22
      8  558  941 (S23E27 L= 72)   530    205        2 ( 465/ 912- 643/ 973)  13.0 E33-E20 S24-S20
      9  347  737 (S39E44 L= 55)   605    129        2 ( 316/ 678- 378/ 795)  53.1 E47-E40 S43-S35
     10 1639  556 (S15W55 L=154)   861    224        2 (1556/ 482-1711/ 636)  71.6 W52-W57 S21-S09
     11 1186  518 (S31W25 L=124)  1133    221        8 (1096/ 464-1266/ 574)  52.1 W21-W28 S36-S26
     12  778  390 (S53W03 L=101)   543    136        1 ( 740/ 356- 837/ 442) -21.8 E01-W09 S54-S49
  

    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