Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1320 1642 (N45W22 L=262)   474     74        1 (1275/1615-1360/1666)   7.1 W18-W26 N42-N47
      2 1587 1603 (N43W49 L=289)  1927    277        3 (1463/1566-1704/1634)  -1.9 W34-W66 N39-N44
      3  558 1374 (N22E34 L=206)   277     53        1 ( 533/1355- 592/1397)  14.0 E35-E31 N20-N23
      4 1601 1133 (N11W38 L=278)  2152    414        1 (1522/ 948-1651/1298) -79.2 W32-W41 N00-N21
      5  275 1197 (N09E55 L=185)   906    114        2 ( 254/1139- 303/1254) -74.1 E57-E51 N05-N12
      6  703 1229 (N14E21 L=218)   258     77        1 ( 683/1187- 720/1254)  90.0 E22-E20 N10-N15
      7  854 1193 (N12E11 L=229)   284    130        1 ( 817/1148- 898/1246) -50.2 E14-E08 N09-N15
      8 1071 1129 (N09W02 L=242)   325    122        4 (1048/1079-1102/1186) -71.6 W00-W04 N06-N12
      9  900  997 (N00E07 L=232)   827    218        1 ( 849/ 905- 947/1081)  90.0 E10-E04 S05-N05
     10  483  861 (S10E35 L=205)   779    241        1 ( 439/ 773- 551/ 961) -57.5 E38-E29 S15-S04
     11 1129  696 (S17W08 L=248)  1834    596        3 (1028/ 478-1226/ 928)  73.3 W02-W12 S32-S02
     12 1681  698 (S16W48 L=288)   878    279        1 (1652/ 611-1722/ 807) -65.2 W44-W54 S21-S08
     13 1770  688 (S16W58 L=298)  1635    220        5 (1726/ 632-1829/ 750) -38.7 W53-W66 S20-S12
     14 1380  286 (S47W39 L=279)   509     86        0 (1340/ 268-1411/ 320)  38.7 W35-W42 S49-S43
     15  813  242 (S55E17 L=223)  1938    208        4 ( 730/ 214- 919/ 277)   9.9 E28-E05 S59-S50
     16 1085  259 (S51W11 L=251)   520     81        1 (1044/ 235-1126/ 273)   7.1 W06-W14 S53-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