Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1264 1713 (N45W38 L=352)  1189    174        1 (1208/1658-1343/1750)   3.6 W30-W46 N41-N46
      2  625 1703 (N60E25 L=289)   711     80        2 ( 583/1660- 653/1738) -63.4 E31-E20 N55-N63
      3 1468 1577 (N31W45 L=359)   443     85        2 (1430/1556-1514/1594) -26.6 W41-W48 N28-N32
      4  993 1449 (N33W06 L=320)  2043    538        3 ( 755/1368-1223/1528) -33.2 E10-W19 N23-N42
      5  255 1403 (N40E63 L=251)   782     87        2 ( 218/1368- 291/1435) -12.5 E73-E55 N38-N42
      6 1599 1334 (N14W44 L=358)   421    132        0 (1552/1276-1630/1371) -71.6 W41-W45 N10-N16
      7  181 1091 (N21E64 L=250)   854    103        2 ( 164/1039- 200/1146)  74.1 E67-E61 N17-N24
      8  275 1074 (N19E53 L=261)   501    134        0 ( 258/1016- 298/1132)  74.1 E54-E50 N15-N22
      9  620 1004 (N11E25 L=289)  1145    376        2 ( 488/ 888- 751/1098) -46.8 E33-E18 N02-N18
     10 1582  904 (S11W33 L=347)   527    135        1 (1564/ 835-1604/ 964)  60.3 W30-W35 S14-S07
     11 1693  942 (S11W42 L=356)   289     94        0 (1655/ 924-1737/ 963) -33.7 W39-W45 S12-S08
     12  916  702 (S11E12 L=302)   799    375        1 ( 857/ 572- 981/ 833)  83.3 E14-E08 S19-S02
     13 1058  681 (S14E04 L=311)   485    211        0 (1021/ 633-1113/ 746) -76.0 E05-W00 S17-S09
     14 1230  472 (S32W04 L=318)  1471    330        4 (1111/ 403-1345/ 590)  27.9 E05-W13 S34-S25
     15  887  465 (S26E19 L=295)   376     79        1 ( 836/ 444- 926/ 481)   0.0 E22-E16 S27-S24
     16  941  340 (S36E20 L=294)   269     48        0 ( 914/ 326- 971/ 355)   0.0 E22-E17 S36-S34
     17  745  224 (S43E44 L=270)  1130    154        2 ( 683/ 187- 822/ 264)  11.3 E54-E34 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