Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1649 1398 (N12W51 L=  7)  1852    350        1 (1531/1310-1779/1531) -62.2 W46-W59 N03-N22
      2  940 1445 (N34W05 L=321)   866    269        2 ( 811/1392-1034/1478) -34.7 E03-W10 N28-N38
      3 1731 1428 (N11W61 L= 16)   510     59        2 (1703/1399-1758/1452)  56.3 W58-W63 N09-N12
      4 1201 1385 (N24W20 L=336)   220     67        0 (1167/1370-1235/1404)   0.0 W17-W22 N22-N24
      5 1299 1252 (N14W22 L=337)   688    280        1 (1221/1177-1375/1322)  13.0 W15-W28 N09-N16
      6  351 1182 (N30E44 L=271)   341     67        0 ( 317/1165- 394/1195) -14.0 E48-E40 N28-N31
      7  694 1089 (N18E19 L=297)  1194    373        0 ( 521/1052- 828/1166)   0.0 E30-E08 N15-N20
      8  997 1039 (N09E01 L=314)   772    149        4 ( 939/1009-1057/1080) -18.4 E04-W02 N06-N11
      9  179  999 (N21E65 L=250)   605     67        2 ( 153/ 972- 205/1021) -15.9 E69-E61 N19-N22
     10 1521  575 (S33W22 L=338)   278     38        1 (1501/ 550-1539/ 597)  63.4 W20-W23 S35-S31
  

    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