Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  611 1653 (N37E41 L=172)   392     50        0 ( 586/1628- 640/1672)  63.4 E43-E39 N35-N39
      2  499 1653 (N36E51 L=162)   642     95        0 ( 459/1632- 544/1670)  31.0 E54-E47 N34-N37
      3  212 1222 (N05E63 L=150)   550     66        0 ( 197/1185- 233/1259)  59.0 E65-E61 N02-N07
      4  663 1133 (N04E24 L=189)   299     58        1 ( 640/1101- 684/1166)  76.0 E24-E22 N02-N06
      5 1098 1109 (N07W04 L=217)   224     72        0 (1062/1080-1137/1137) -31.0 W01-W06 N05-N08
      6 1104  890 (S06W06 L=219)   847    257        1 (1038/ 767-1147/ 996) -58.4 W01-W10 S13-N00
      7  766  872 (S11E15 L=199)   299     86        2 ( 724/ 858- 809/ 893)  18.4 E17-E11 S11-S08
      8 1193  682 (S18W14 L=227)   259     49        0 (1169/ 660-1223/ 699)  18.4 W12-W16 S19-S16
  

    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