Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1354 1571 (N35W05 L=114)  1928    418        1 (1220/1372-1451/1706) -34.8 E08-W15 N24-N41
      2  238 1433 (S01E69 L= 40)  4391    552        3 ( 125/1197- 335/1599)  83.2 E77-E65 S15-N11
      3  615 1555 (N13E40 L= 69)   537    145        2 ( 553/1528- 687/1595)  45.0 E44-E36 N10-N17
      4  935 1465 (N15E17 L= 92)   365    102        0 ( 883/1435- 987/1485)   0.0 E20-E13 N14-N16
      5 1023 1280 (N07E07 L=102)   223     78        1 ( 998/1247-1049/1316) -36.9 E09-E04 N05-N08
      6  292 1162 (S16E51 L= 58)  1025    268        2 ( 249/1046- 376/1268)  71.6 E54-E46 S23-S08
      7  173 1163 (S18E64 L= 45)  1203    143        3 ( 139/1107- 196/1216)  81.9 E67-E60 S21-S14
      8 1339  594 (S22W32 L=141)   594    155        2 (1262/ 577-1416/ 616)  29.1 W26-W36 S23-S18
      9 1526  391 (S26W55 L=164)   652    112        0 (1476/ 367-1587/ 413)  40.6 W52-W59 S28-S22
  

    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