Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1049 1283 (N19W09 L=325)   628    128        0 ( 990/1255-1123/1312) -39.8 W05-W12 N16-N22
      2  165 1110 (N31E59 L=257)  2148    321        1 ( 109/1021- 218/1196) -53.1 E65-E54 N24-N36
      3 1861 1052 (S18W57 L= 13)   414     72        1 (1845/1024-1879/1081)  21.8 W54-W59 S19-S16
      4  816  953 (N07E13 L=303)   872    234        1 ( 702/ 914- 939/ 981) -24.8 E19-E06 N03-N10
      5  456  918 (N14E37 L=279)   671    180        0 ( 391/ 861- 522/ 959)   9.5 E43-E31 N11-N15
      6 1284  872 (S10W10 L=327)   381    101        0 (1232/ 841-1321/ 914)   8.1 W06-W13 S10-S08
      7  304  853 (N13E51 L=265)   389    109        0 ( 259/ 833- 332/ 871) -18.4 E54-E48 N11-N14
      8 1491  825 (S18W22 L=338)   241     44        0 (1470/ 791-1510/ 851)  45.0 W20-W23 S20-S17
      9  186  781 (N11E66 L=250)   470     79        0 ( 167/ 750- 212/ 822) -68.2 E67-E63 N08-N13
     10  988  491 (S24E18 L=299)  1558    352        2 ( 944/ 346-1043/ 636)  45.0 E26-E09 S32-S16
     11  365  450 (S12E64 L=252)  7037    814       12 ( 283/ 351- 466/ 542) -50.2 E75-E56 S19-S05
     12  481  445 (S14E53 L=263)   537     74        0 ( 454/ 422- 510/ 472)  26.6 E55-E49 S15-S12
  

    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