Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  968 1521 (N24E17 L=285)  1370    481        0 ( 834/1344-1094/1681) -28.2 E32-E04 N16-N31
      2 1260 1594 (N37E01 L=301)   992    204        2 (1209/1499-1309/1664) -39.3 E07-W04 N31-N41
      3 1644 1617 (N55W37 L=339)   679     58        0 (1622/1585-1676/1637)  32.5 W32-W44 N51-N58
      4  404  889 (S27E38 L=265)  1290    348        0 ( 308/ 771- 533/ 954)  39.8 E46-E27 S36-S20
      5  197  853 (S34E61 L=241)   703     97        3 ( 176/ 811- 226/ 893)  59.0 E64-E57 S36-S31
      6  894  616 (S31W03 L=306)  1079    277        0 ( 749/ 595-1013/ 642)  26.6 E06-W10 S34-S26
      7 1138  563 (S27W20 L=323)   299     60        0 (1117/ 527-1162/ 591) -76.0 W19-W21 S29-S25
      8 1444  505 (S22W43 L=345)   285     28        1 (1422/ 489-1470/ 523)   0.0 W40-W45 S22-S20
      9 1353  467 (S26W39 L=341)   819    112        1 (1318/ 412-1383/ 520) -69.4 W35-W41 S30-S22
     10  653  325 (S61E07 L=295)  1377    245        2 ( 556/ 257- 747/ 405) -12.1 E20-W08 S64-S56
     11 1439  348 (S31W54 L=357)   962    234        1 (1357/ 302-1524/ 384)  13.2 W47-W64 S33-S27
     12  911  216 (S60W29 L=331)   586     62        0 ( 880/ 189- 946/ 236)   0.0 W23-W35 S61-S58
  

    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