Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1147 1874 (N60W51 L=286)   485     59        0 (1128/1847-1176/1904)   9.5 W44-W62 N57-N61
      2 1068 1610 (N42W18 L=253)  1378    270        4 ( 934/1576-1188/1635) -14.0 W07-W27 N38-N44
      3  474 1494 (N46E36 L=199)  1977    374        3 ( 334/1427- 639/1567)   8.9 E52-E19 N43-N48
      4 1826 1265 (N03W62 L=297)   575     60        2 (1795/1237-1849/1289) -31.0 W59-W65 N00-N04
      5  229  941 (N13E58 L=177)   661     64        0 ( 204/ 910- 250/ 989)  51.3 E61-E56 N11-N16
      6  505  939 (N10E34 L=201)   279     88        0 ( 481/ 901- 527/ 979) -78.7 E35-E32 N07-N12
      7 1156  714 (S14W02 L=237)   283     89        0 (1139/ 669-1176/ 758)  68.2 W00-W03 S16-S11
      8 1222  594 (S23W04 L=239)   520    160        0 (1169/ 559-1299/ 646) -45.0 W01-W08 S26-S18
      9  575  499 (S18E41 L=194)   853    220        2 ( 498/ 437- 616/ 585)  29.7 E49-E35 S21-S12
     10 1456  397 (S44W19 L=254)  3078    537        4 (1231/ 285-1674/ 520)   8.1 E02-W41 S47-S39
     11  857  337 (S34E28 L=207)   319     84        0 ( 826/ 307- 880/ 378)  29.7 E31-E24 S35-S30
  

    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