Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1250 1812 (N55E14 L=243)   466     56        0 (1215/1801-1286/1824)  33.7 E17-E10 N52-N56
      2  891 1787 (N38E37 L=220)   818    165        0 ( 806/1766- 956/1803)  18.4 E44-E31 N35-N40
      3  299 1524 (N07E69 L=189)  2389    172        5 ( 243/1463- 368/1576)  48.4 E76-E63 N02-N11
      4 1335 1342 (N22W09 L=267)   450    155        0 (1298/1304-1387/1397) -23.2 W06-W13 N20-N24
      5  229 1291 (S08E63 L=194)   629     83        1 ( 186/1271- 262/1316)  26.6 E68-E58 S10-S06
      6 1863  884 (N12W63 L=320)  1070    132        1 (1832/ 832-1890/ 943)  69.4 W59-W66 N08-N16
      7 1351  570 (S21W34 L=291)   300     98        1 (1334/ 517-1363/ 615) -56.3 W31-W36 S23-S17
      8  759  547 (S40E03 L=254)   424    118        0 ( 732/ 505- 807/ 598) -33.7 E05-W02 S42-S37
      9 1403  501 (S23W41 L=298)   314     49        1 (1377/ 486-1433/ 520)   0.0 W38-W42 S24-S21
  

    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