Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1306 1647 (N45W21 L=262)   555     75        0 (1267/1619-1348/1671)  20.6 W17-W25 N42-N47
      2 1576 1602 (N43W47 L=289)  2017    264        2 (1444/1560-1691/1632)   0.0 W32-W64 N39-N44
      3  240 1446 (N24E71 L=171)   869     79        0 ( 214/1410- 257/1488) -51.3 E74-E66 N21-N26
      4  541 1376 (N22E35 L=206)   277     60        1 ( 514/1358- 574/1399)  11.3 E37-E32 N20-N23
      5 1582 1134 (N11W36 L=278)  2424    416        4 (1505/ 951-1630/1309) -74.1 W31-W40 N00-N21
      6  684 1237 (N14E23 L=219)   232     34        1 ( 661/1214- 700/1261) -33.7 E24-E21 N12-N15
      7  831 1199 (N12E13 L=228)   308    119        0 ( 798/1155- 875/1249) -39.8 E15-E09 N09-N15
      8  280 1143 (N06E53 L=188)  1341    192        2 ( 253/1042- 315/1226) -63.4 E57-E49 N00-N10
      9 1049 1113 (N08W01 L=242)   272     75        0 (1027/1081-1078/1148) -63.4 W00-W03 N06-N10
     10  875  995 (N00E09 L=232)   945    244        3 ( 811/ 888- 920/1084)  79.7 E12-E06 S06-N05
     11 1119  727 (S15W07 L=249)  2133    658        2 (1011/ 486-1247/ 993)  70.0 W01-W13 S32-N01
     12  465  862 (S10E36 L=205)   916    246        1 ( 417/ 774- 536/ 966) -54.0 E39-E31 S15-S04
     13 1659  776 (S11W45 L=286)   326     75        0 (1643/ 737-1684/ 811)  90.0 W43-W46 S13-S08
     14 1753  692 (S16W56 L=297)  1575    183        4 (1712/ 627-1801/ 761) -56.3 W51-W61 S20-S11
     15 1673  645 (S19W49 L=290)   553    126        0 (1636/ 610-1722/ 694) -21.8 W44-W54 S21-S15
     16 1360  282 (S48W38 L=279)   446     83        0 (1328/ 264-1393/ 318)  45.0 W34-W40 S49-S44
     17  866  248 (S54E11 L=230)  2694    412        6 ( 717/ 218-1088/ 281)   8.5 E29-W11 S58-S49
  

    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