Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  274 1441 (N31E64 L=324)  2901    226        7 ( 213/1353- 322/1522) -23.7 E82-E53 N25-N36
      2 1002 1473 (N33W00 L= 29)   284     86        0 ( 981/1437-1033/1510) -68.2 E01-W02 N30-N35
      3  890 1420 (N30E08 L= 21)   319    124        0 ( 843/1375- 924/1446) -45.0 E11-E05 N26-N31
      4  796 1308 (N22E14 L= 15)   203     46        1 ( 771/1287- 820/1325) -33.7 E15-E12 N21-N23
      5  531 1270 (N20E33 L=356)  1014    231        2 ( 430/1236- 626/1302)  -3.6 E41-E25 N18-N22
      6 1132 1249 (N18W08 L= 36)   199     65        0 (1094/1237-1159/1267)  14.0 W05-W09 N16-N18
      7  883 1163 (N13E08 L= 20)   216     77        1 ( 843/1148- 913/1184)  14.0 E10-E06 N12-N14
      8 1576 1068 (N04W36 L= 65)   302     80        0 (1542/1040-1608/1101)  21.8 W33-W38 N02-N06
      9 1638  973 (S02W40 L= 69)   745    230        0 (1615/ 867-1658/1078) -85.2 W38-W41 S08-N04
     10 1280  971 (S00W16 L= 44)   274     78        0 (1264/ 935-1294/1020)  76.0 W14-W16 S02-N02
     11  823  810 (S08E13 L= 15)   340    115        0 ( 775/ 769- 874/ 841)  29.7 E16-E09 S10-S06
     12 1554  805 (S12W34 L= 63)   671    196        1 (1478/ 746-1630/ 839)  24.4 W28-W39 S15-S10
  

    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