Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  774 1541 (N26E20 L=198)   890    236        1 ( 688/1461- 831/1604) -33.7 E27-E15 N20-N29
      2  313 1410 (N16E53 L=166)  3613    435       16 ( 191/1303- 437/1519)  23.2 E67-E44 N08-N22
      3  434 1442 (N18E43 L=176)   315     71        1 ( 419/1412- 455/1476) -78.7 E44-E41 N15-N20
      4  189 1300 (N09E64 L=155)  2275    178        2 ( 137/1217- 247/1374) -35.5 E74-E57 N03-N13
      5 1302 1317 (N15W15 L=234)   219     38        0 (1286/1292-1320/1342) -90.0 W14-W16 N13-N16
      6  664 1085 (S03E22 L=197)   431    139        1 ( 611/1029- 699/1144)  50.2 E25-E20 S06-N00
      7  125 1079 (S04E69 L=150)   493     71        0 (  93/1060- 158/1095)  11.3 E74-E63 S05-S03
      8  487 1005 (S09E34 L=185)   430    129        0 ( 435/ 958- 513/1047) -39.8 E37-E31 S11-S06
      9  370  979 (S10E43 L=176)   320     60        0 ( 351/ 949- 394/1005)  53.1 E44-E40 S12-S08
     10  960  963 (S09E03 L=215)   289    108        0 ( 910/ 942-1012/ 981)  -9.5 E06-W00 S09-S07
     11  399  880 (S16E41 L=178)  1119    260        3 ( 338/ 819- 454/ 953)  48.8 E45-E36 S19-S11
     12  248  869 (S17E55 L=164)  2617    422        6 ( 170/ 765- 317/ 949)  32.0 E66-E47 S23-S12
     13 1059  815 (S17W03 L=222)   255    111        0 (1008/ 799-1114/ 839)   9.5 W00-W06 S17-S15
     14 1196  528 (S34W16 L=235)   702    182        0 (1105/ 477-1261/ 563) -18.4 W09-W21 S36-S32
     15 1440  499 (S33W35 L=254)   305     76        2 (1391/ 482-1482/ 518)  23.2 W31-W38 S34-S31
     16  874  472 (S40E07 L=212)   674    175        1 ( 802/ 443- 968/ 505)  14.0 E12-W00 S42-S37
     17 1354  256 (S52W44 L=263)  1776    210        3 (1246/ 223-1456/ 289)   4.2 W31-W59 S54-S48
  

    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