Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  994 1878 (N47E37 L=111)   582     64        1 ( 980/1843-1018/1908) -33.7 E41-E32 N43-N49
      2  779 1850 (N38E49 L= 99)   870     79        1 ( 733/1830- 816/1874)   7.1 E52-E44 N35-N40
      3 1199 1628 (N33E07 L=141)   674     92        2 (1163/1580-1239/1666) -23.2 E10-E03 N30-N35
      4  632 1473 (N09E35 L=113)   601     97        2 ( 589/1444- 689/1498)  45.0 E38-E31 N06-N10
      5  434 1440 (N03E48 L=100)   484    102        1 ( 401/1403- 469/1471)  53.1 E49-E45 N00-N05
      6  987 1270 (N05E08 L=139)   480    170        0 ( 924/1230-1033/1335)  69.4 E11-E06 N02-N10
      7 1769 1261 (N28W43 L=191)   514    116        0 (1726/1206-1796/1302) -29.7 W38-W46 N25-N30
      8  557 1227 (S07E32 L=116)   422    109        1 ( 517/1180- 589/1272) -26.6 E36-E29 S08-S05
      9  168 1224 (S14E67 L= 80)   542     64        2 ( 144/1195- 199/1244)  31.0 E71-E63 S15-S12
     10  643 1060 (S14E23 L=125)   814    282        2 ( 561/1025- 740/1095)  24.4 E27-E16 S17-S11
     11  788  864 (S22E09 L=138)   190     62        0 ( 757/ 853- 825/ 878)  26.6 E11-E07 S22-S20
     12 1492  648 (S14W39 L=187)  1392    304        2 (1352/ 604-1625/ 710)  40.9 W31-W46 S20-S07
  

    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