Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  970 1522 (N37E00 L=297)   824    275        0 ( 838/1464-1079/1577)  15.5 E10-W08 N33-N40
      2 1774 1449 (N24W66 L=  3)  1012    124        3 (1736/1407-1814/1489) -50.2 W61-W72 N20-N26
      3 1842 1424 (N21W76 L= 14)  1494     71        1 (1821/1395-1864/1463)  76.0 W71-W83 N19-N23
      4  361 1214 (N19E46 L=251)   500     92        0 ( 327/1176- 390/1254) -39.8 E49-E43 N16-N21
      5 1181 1164 (N12W11 L=308)   411    154        0 (1116/1131-1247/1186) -14.0 W06-W14 N09-N13
      6  194 1108 (N12E63 L=234)   555     82        0 ( 171/1072- 221/1151) -33.7 E66-E59 N10-N14
      7  867  894 (S02E10 L=287)  1744    489        1 ( 789/ 683- 993/1104) -83.2 E14-E03 S15-N10
      8  536 1055 (N09E31 L=266)   271     78        0 ( 519/1019- 567/1090)  68.2 E32-E28 N06-N11
      9 1641 1023 (N00W41 L=338)   267     88        0 (1623/ 985-1656/1074)  53.1 W39-W42 S02-N03
     10 1678  969 (S04W43 L=340)   756    212        1 (1602/ 926-1742/1050)  18.4 W37-W49 S05-N00
     11 1698  832 (S12W45 L=342)   481    115        0 (1675/ 773-1733/ 871)  59.0 W43-W48 S15-S10
     12  402  669 (S15E46 L=251)   612    113        1 ( 365/ 623- 435/ 722) -63.4 E48-E43 S17-S11
     13 1479  593 (S26W29 L=326)   439     88        1 (1433/ 565-1522/ 628)  23.2 W25-W32 S27-S23
     14  403  595 (S20E48 L=249)   425     69        1 ( 374/ 572- 427/ 623) -33.7 E50-E45 S21-S17
     15  679  466 (S29E29 L=268)   258     61        1 ( 667/ 427- 694/ 499)  63.4 E30-E27 S32-S26
     16 1253  345 (S43W14 L=311)  1639    181        7 (1206/ 278-1305/ 405)  63.4 W09-W18 S49-S37
     17  905  287 (S45E17 L=280)   431     61        0 ( 877/ 260- 927/ 319)  45.0 E20-E14 S47-S42
     18  706  252 (S48E38 L=259)  2808    287        6 ( 573/ 210- 843/ 292)   4.6 E63-E22 S52-S44
  

    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