Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  613 1761 (N68E33 L=194)   564     57        1 ( 586/1742- 649/1781)  -4.4 E39-E25 N65-N69
      2 1013 1639 (N46W14 L=241)   255     29        0 ( 991/1623-1034/1657)  11.3 W11-W16 N44-N47
      3 1631 1420 (N17W51 L=278)  1740    397        0 (1468/1351-1816/1472) -26.6 W39-W68 N09-N24
      4 1212 1414 (N27W20 L=247)   433    147        2 (1161/1375-1246/1452) -59.0 W17-W22 N24-N29
      5  541 1180 (N23E29 L=198)  1049    217        0 ( 500/1075- 576/1279) -81.9 E32-E26 N16-N30
      6 1637 1164 (N03W42 L=269)   733    170        2 (1602/1085-1659/1233)  76.0 W39-W44 S01-N07
      7 1359 1163 (N09W23 L=250)   475    188        1 (1286/1128-1427/1190)  -5.7 W17-W27 N05-N10
      8  471  660 (S07E43 L=184)   564    173        1 ( 400/ 610- 538/ 710) -45.0 E47-E38 S11-S03
      9  667  632 (S11E29 L=198)   228     96        0 ( 653/ 603- 700/ 672) -76.0 E30-E27 S13-S08
     10  912  393 (S32E21 L=206)  1243    355        1 ( 792/ 311-1027/ 508)  22.6 E33-E09 S37-S24
  

    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