Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1589 1680 (N25W67 L= 63)   517     80        1 (1540/1671-1628/1694) -16.7 W62-W72 N22-N26
      2 1526 1606 (N23W54 L= 50)   739     97        3 (1484/1563-1569/1640)  -5.7 W49-W59 N21-N25
      3 1388 1563 (N26W41 L= 37)   568    120        0 (1352/1506-1416/1609)  60.3 W38-W43 N21-N28
      4 1354 1419 (N19W32 L= 27)   278     57        0 (1331/1394-1377/1443)  14.0 W29-W33 N18-N20
      5 1134 1399 (N24W17 L= 13)   276     72        0 (1095/1382-1178/1418) -45.0 W15-W19 N22-N26
      6  990 1296 (N23W06 L=  1)   731    245        2 ( 871/1238-1070/1332) -41.6 E01-W08 N17-N26
      7  710 1203 (N25E14 L=342)   356    106        1 ( 680/1155- 736/1251) -90.0 E15-E12 N21-N27
      8  274 1204 (N35E51 L=304)   809    148        1 ( 222/1164- 340/1249) -21.8 E59-E43 N32-N38
      9 1909 1074 (S18W67 L= 62)  1409    134        2 (1848/1034-1956/1107) -21.8 W57-W78 S22-S14
     10 1635 1027 (S11W37 L= 32)   617    148        0 (1595/ 968-1690/1083)  18.4 W32-W42 S13-S08
     11 1249  774 (S13W06 L=  2)   924    306        0 (1189/ 655-1308/ 910)  42.5 W00-W12 S18-S06
     12  500  825 (N08E37 L=319)   733    166        1 ( 448/ 771- 558/ 894) -66.0 E39-E33 N03-N12
     13 1567  811 (S21W27 L= 23)   274     64        0 (1547/ 778-1584/ 842) -90.0 W25-W28 S23-S18
     14 1264  520 (S29W00 L=356)   474     60        2 (1239/ 485-1295/ 554)  63.4 E02-W02 S31-S27
     15  985  524 (S21E16 L=339)   218     73        0 ( 960/ 506-1019/ 550) -45.0 E17-E14 S22-S18
     16  813  420 (S23E31 L=325)   537    189        0 ( 769/ 348- 849/ 484)  36.0 E38-E27 S27-S19
     17  594  275 (S28E59 L=297)   818    121        1 ( 561/ 231- 620/ 330)  23.6 E68-E52 S31-S24
     18  802  229 (S37E46 L=310)   520     85        0 ( 776/ 190- 830/ 270)  58.0 E49-E42 S41-S33
  

    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