Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1389 1361 (N16W31 L=211)   413    122        1 (1342/1318-1425/1396) -71.6 W29-W33 N13-N19
      2  980 1359 (N27W06 L=186)   371    151        1 ( 925/1310-1045/1382)   0.0 W01-W10 N25-N28
      3  536 1309 (N35E26 L=154)  1297    376        1 ( 404/1254- 715/1352) -13.5 E37-E12 N30-N38
      4 1794 1266 (S01W58 L=238)   653     65        2 (1758/1233-1830/1288) -14.0 W54-W62 S03-N00
      5  350 1237 (N35E44 L=136)   455     79        0 ( 328/1197- 376/1274)  45.0 E45-E40 N32-N36
      6 1434 1104 (N01W26 L=206)   927    302        2 (1389/ 976-1475/1238)  61.7 W23-W30 S06-N07
      7 1874  980 (S20W59 L=239)  1093    187        0 (1824/ 896-1926/1040)  18.4 W52-W69 S23-S17
      8 1346  660 (S22W10 L=190)   858    303        0 (1228/ 586-1424/ 763)  17.4 W00-W16 S24-S18
      9  492  474 (S14E51 L=129)   429     52        2 ( 463/ 448- 521/ 501) -63.4 E52-E48 S15-S11
     10  930  399 (S28E24 L=156)   279     81        0 ( 892/ 382- 973/ 410)  -9.5 E26-E20 S28-S26
  

    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