Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1112 1705 (N42W31 L=130)   341    121        0 (1064/1681-1178/1730) -12.5 W26-W36 N40-N42
      2 1539 1448 (N14W44 L=144)  1587    455        3 (1368/1372-1728/1512) -39.0 W35-W56 N04-N21
      3 1048 1459 (N29W14 L=114)   666    261        1 ( 914/1432-1169/1483) -26.6 W06-W22 N25-N34
      4 1843 1384 (N00W70 L=170)   905    255        0 (1808/1359-1886/1401)  -7.6 W64-W79 S01-N01
      5 1359 1357 (N14W29 L=128)   380    103        2 (1315/1335-1411/1381) -38.7 W26-W31 N12-N16
      6 1847 1123 (S13W57 L=157)  2750    400        6 (1804/ 946-1894/1263)  75.3 W53-W64 S24-S04
      7 1139 1129 (N08W09 L=109)   187    102        1 (1101/1094-1165/1158)  11.3 W06-W11 N06-N09
      8 1177 1045 (N02W09 L=108)   259    100        1 (1129/1029-1234/1065) -38.7 W06-W11 N00-N04
      9  323  955 (N18E47 L= 53)   854    211        4 ( 271/ 907- 407/1002) -28.6 E50-E39 N14-N21
     10 1651  923 (S18W35 L=135)   908    256        3 (1568/ 835-1708/ 989) -69.0 W30-W38 S24-S11
     11 1163  652 (S19E02 L= 98)   558    215       -1 (1120/ 609-1215/ 689)   9.5 E05-W01 S21-S16
     12 1541  629 (S33W22 L=121)   275     80        1 (1518/ 597-1565/ 671)  53.1 W19-W23 S34-S29
     13 1293  499 (S33W01 L=101)   876    234        1 (1271/ 424-1322/ 620)  74.1 E01-W03 S38-S24
     14  745  438 (S21E34 L= 65)   520    141        0 ( 718/ 383- 790/ 477)  35.5 E37-E30 S24-S19
     15 1483  394 (S48W13 L=113)  2220    352        8 (1409/ 342-1562/ 459)   8.1 W06-W21 S54-S42
     16  500  319 (S24E59 L= 40)  2077    294        7 ( 471/ 258- 540/ 384)  48.0 E66-E54 S28-S18
     17  770  340 (S29E38 L= 62)   491     90        0 ( 736/ 310- 812/ 382) -80.5 E39-E36 S31-S25
  

    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