Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1510 1635 (N43W43 L=131)  2114    361        7 (1447/1568-1555/1707)  32.0 W35-W52 N38-N49
      2 1065 1499 (N34W03 L= 90)   503    154        1 (1045/1434-1087/1562)  90.0 W01-W04 N29-N38
      3  716 1509 (N35E23 L= 64)   239     78        0 ( 702/1472- 730/1544) -68.2 E25-E21 N32-N37
      4 1666 1320 (N21W46 L=134)   472    126        2 (1617/1275-1709/1345) -33.7 W42-W49 N17-N22
      5 1038 1168 (N12W01 L= 88)   309    189        0 (1012/1088-1064/1266)  79.7 W00-W02 N07-N18
      6 1256 1188 (N14W15 L=102)   538    146        1 (1213/1119-1284/1262) -66.0 W12-W16 N09-N18
      7 1667 1232 (N15W45 L=132)   785    199        1 (1608/1190-1733/1262)   9.5 W39-W51 N12-N17
      8 1691 1110 (N08W46 L=133)   825    249        3 (1610/1078-1782/1167) -18.4 W39-W54 N05-N11
      9 1446 1116 (N09W27 L=114)   264    106        2 (1407/1067-1489/1153) -39.8 W24-W30 N05-N11
     10  644  983 (N01E24 L= 63)   353    112        1 ( 592/ 956- 700/1019)  15.9 E27-E20 N00-N03
     11 1714  704 (S18W50 L=137)  1468    327        2 (1637/ 576-1757/ 862)  75.3 W45-W53 S26-S07
     12  269  813 (S11E55 L= 32)   376     82        1 ( 241/ 798- 299/ 839)   0.0 E58-E52 S12-S09
     13  984  736 (S14E03 L= 85)   281    117        1 ( 920/ 716-1039/ 757)   0.0 E06-W00 S15-S13
     14  608  585 (S25E29 L= 58)   390     98        0 ( 559/ 569- 652/ 606)   7.1 E33-E25 S26-S23
     15 1250  417 (S37W18 L=105)   500    111        2 (1221/ 379-1278/ 460) -66.8 W15-W20 S40-S33
     16 1158  348 (S43W11 L= 99)  1251    174        5 (1102/ 298-1227/ 411)   5.7 W06-W16 S47-S37
     17  846  364 (S42E15 L= 73)   512     73        1 ( 818/ 341- 884/ 400) -21.8 E17-E11 S43-S38
     18  713  347 (S43E27 L= 60)   546    146        2 ( 671/ 312- 751/ 383)  31.0 E32-E22 S46-S40
  

    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