Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1367 1772 (N39W53 L=198)  1704    294        2 (1242/1707-1469/1836)   0.0 W37-W74 N37-N40
      2  680 1416 (N34E13 L=132)   412    118        1 ( 629/1397- 749/1434) -20.6 E16-E08 N31-N35
      3 1774 1340 (N03W59 L=204)   493     79        0 (1744/1318-1803/1361) -26.6 W55-W61 N01-N04
      4  913 1320 (N22W00 L=146)   203     71        1 ( 876/1306- 948/1334)   0.0 E02-W02 N21-N23
      5 1402 1263 (N08W28 L=173)   267     75        0 (1375/1231-1426/1292) -63.4 W26-W28 N05-N09
      6  127 1151 (N29E71 L= 74)  1114    112        1 ( 101/1102- 158/1202) -25.0 E81-E64 N25-N32
      7  969 1178 (N13W00 L=146)   189    114        0 ( 918/1161-1020/1202) -11.3 E02-W03 N12-N13
      8 1821 1181 (S07W58 L=203)   440     44        0 (1794/1160-1842/1202)  -9.5 W54-W60 S08-S05
      9 1865 1081 (S14W60 L=206)  2616    193        6 (1808/1010-1916/1175)  23.6 W52-W69 S18-S08
     10  413 1057 (N18E38 L=108)  1284    349        0 ( 269/1013- 537/1145)   2.5 E50-E27 N15-N20
     11 1639 1038 (S10W38 L=184)   270    105        0 (1627/ 992-1655/1091)  68.2 W36-W39 S12-S07
     12 1736  877 (S22W44 L=189)   401     78        0 (1713/ 838-1754/ 919)  51.3 W41-W45 S24-S19
     13 1589  807 (S23W31 L=176)   835    206        1 (1518/ 742-1637/ 893)  35.0 W24-W35 S25-S18
     14 1117  755 (S14E01 L=145)   354    211        0 (1055/ 676-1179/ 835)  35.0 E06-W04 S17-S10
     15 1393  731 (S22W16 L=161)   600    197        0 (1358/ 635-1424/ 818)  51.3 W11-W19 S27-S17
     16 1637  715 (S30W34 L=179)   336     72        1 (1607/ 685-1666/ 753)  36.9 W31-W36 S31-S27
     17 1036  552 (S25E11 L=135)   204     76        0 (1012/ 528-1060/ 592)  36.9 E12-E08 S25-S22
     18  726  367 (S30E38 L=107)  1054    328        0 ( 705/ 219- 751/ 520)  43.7 E52-E29 S40-S19
  

    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