Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  652 1598 (N18E39 L=184)   471    127        1 ( 621/1554- 684/1663) -26.6 E44-E34 N15-N22
      2  358 1507 (N07E58 L=165)  1326    233        5 ( 320/1407- 387/1593) -69.0 E64-E52 N00-N13
      3  544 1478 (N09E41 L=182)   590    158        1 ( 511/1408- 592/1547) -37.9 E45-E35 N05-N12
      4  432 1247 (S07E42 L=181)   264     38        1 ( 410/1223- 456/1269) -33.7 E43-E39 S08-S05
      5  342 1193 (S11E48 L=175)   559    142        0 ( 298/1157- 398/1219)   6.3 E51-E42 S14-S09
      6  954 1140 (S02E07 L=216)   198     37        0 ( 931/1122- 977/1156)   0.0 E08-E04 S03-S01
      7  513  938 (S23E30 L=193)   274     81        0 ( 489/ 906- 539/ 977) -31.0 E32-E27 S24-S21
      8  244  870 (S32E55 L=168)   400     53        0 ( 227/ 841- 267/ 893) -71.6 E56-E51 S33-S30
      9  228  785 (S38E58 L=164)   709     79        0 ( 186/ 767- 261/ 812)   9.5 E65-E53 S38-S35
  

    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