Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  532 1124 (N24E29 L=318)   560    197        0 ( 456/1093- 629/1150) -26.6 E34-E22 N20-N26
      2 1541  984 (S08W30 L= 16)   482    130        0 (1496/ 936-1600/1007)   6.3 W25-W34 S09-S06
      3  842  867 (N02E14 L=333)   829    131        3 ( 815/ 750- 876/ 968)  71.6 E17-E11 S04-N08
      4 1680  806 (S23W37 L= 23)   864    244        0 (1632/ 688-1709/ 933)  60.3 W31-W40 S29-S15
      5 1379  828 (S13W16 L=  3)   658    173        1 (1320/ 748-1409/ 913)  45.0 W11-W19 S16-S08
      6 1392  618 (S26W13 L=359)   734    194        1 (1304/ 549-1475/ 711)  19.7 W05-W19 S28-S22
      7 1712  661 (S34W38 L= 25)   606    113        0 (1664/ 628-1755/ 699)   6.3 W33-W42 S35-S32
      8 1004  484 (S24E15 L=332)   309     80        0 ( 973/ 453-1038/ 518)  18.4 E18-E12 S25-S22
      9  957  275 (S39E28 L=319)   519    172        1 ( 935/ 221- 989/ 335)  42.0 E32-E22 S43-S34
  

    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