Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1033 1475 (N23E02 L=218)   351     68        0 ( 987/1452-1069/1497) -11.3 E05-W00 N21-N24
      2 1862 1327 (N21W61 L=282)   473     40        1 (1845/1302-1877/1347)  56.3 W59-W63 N19-N22
      3 1471 1193 (N08W27 L=247)   613    214        0 (1393/1110-1563/1277) -33.7 W21-W33 N04-N13
      4  136 1192 (N04E69 L=151)   536    114        4 ( 117/1168- 161/1220)  56.3 E72-E65 N02-N05
      5 1759 1183 (N10W49 L=270)   398    128        2 (1711/1168-1801/1204)   6.3 W44-W53 N09-N12
      6  650 1058 (S04E23 L=197)   795    206        4 ( 628/ 964- 678/1169) -84.8 E24-E21 S09-N02
      7 1901 1086 (N06W65 L=285)  1009    234        0 (1868/1047-1932/1135)  50.2 W60-W69 N03-N09
      8 1387  747 (S19W25 L=245)  1698    361        0 (1232/ 633-1573/ 813)  27.6 W16-W37 S26-S13
      9  247  649 (S30E58 L=163)   613    144        1 ( 231/ 600- 273/ 688)  39.8 E60-E54 S32-S27
     10 1149  539 (S34W12 L=233)   268     91        0 (1096/ 522-1193/ 565)  33.7 W08-W14 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