Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  700 1729 (N31E44 L=352)  1246    240        0 ( 656/1672- 760/1793) -25.0 E51-E36 N28-N35
      2  646 1589 (N20E40 L=357)  1543    322        4 ( 534/1510- 714/1634)  -6.3 E50-E32 N16-N22
      3 1191 1484 (N28E02 L= 35)  1388    433        1 (1021/1348-1386/1619) -14.9 E17-W13 N23-N33
      4  466 1513 (N12E50 L=347)   706    290        2 ( 442/1387- 515/1604) -70.3 E53-E44 N04-N18
      5 1605 1436 (N37W29 L= 66)   535    125        0 (1551/1394-1662/1470)  -5.7 W24-W34 N35-N38
      6  813 1374 (N11E21 L= 15)   316     80        0 ( 786/1330- 841/1415) -38.7 E23-E18 N09-N13
      7 1748 1293 (N32W43 L= 80)   543    102        0 (1705/1254-1781/1318)  -8.1 W38-W46 N30-N33
      8  738 1176 (S02E20 L= 16)   433    117        0 ( 672/1157- 792/1208)   7.1 E24-E16 S02-N00
      9  699  935 (S16E17 L= 20)   231     90        1 ( 671/ 887- 730/ 986) -45.0 E19-E14 S18-S13
     10  859  502 (S39W04 L= 41)   352     61        1 ( 823/ 482- 896/ 529)   9.5 W01-W07 S39-S37
  

    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