Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1356 1735 (N47W44 L=131)  1726    289        3 (1309/1643-1400/1815)  49.8 W38-W52 N39-N52
      2  937 1567 (N42W02 L= 90)   265     55        1 ( 921/1543- 958/1595) -90.0 W00-W03 N40-N44
      3 1125 1502 (N34W15 L=102)   958    196        1 (1047/1422-1199/1570) -65.6 W09-W17 N28-N39
      4  880 1443 (N34E04 L= 83)   345     87        0 ( 861/1405- 900/1487) -90.0 E05-E03 N31-N36
      5 1707 1253 (N09W50 L=138)  1348    353        4 (1581/1167-1817/1336) -37.4 W42-W60 N02-N16
      6 1281 1203 (N13W18 L=106)   288    108        1 (1222/1172-1327/1230)   8.1 W14-W21 N12-N14
      7  364 1042 (N13E45 L= 43)   671    156        1 ( 309/1000- 415/1100) -37.9 E49-E40 N10-N17
      8  413  706 (S07E46 L= 42)   351     52        1 ( 390/ 683- 436/ 729)  26.6 E47-E43 S08-S06
      9 1562  579 (S31W32 L=120)  2325    416        5 (1390/ 491-1713/ 708)  19.1 W18-W44 S34-S24
     10  787  491 (S25E23 L= 64)   854    181        3 ( 717/ 454- 854/ 546)  15.3 E29-E18 S27-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