Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1010 1882 (N58E37 L= 37)  1102    112        2 ( 952/1851-1059/1908)  58.0 E43-E32 N52-N62
      2  708 1557 (N24E34 L= 41)   447     97        1 ( 661/1534- 757/1584)   0.0 E37-E29 N22-N25
      3  600 1459 (N15E37 L= 37)   335     97        0 ( 587/1414- 622/1497) -78.7 E38-E35 N12-N17
      4 1336 1209 (N16W15 L= 89)  1203    481        4 (1186/1065-1461/1303) -20.9 W03-W25 N09-N19
      5 1529 1228 (N21W28 L=102)   228     40        1 (1508/1206-1549/1250)  56.3 W26-W29 N19-N22
      6 1555 1163 (N18W31 L=105)   414     98        1 (1508/1129-1602/1183)  23.2 W27-W34 N14-N19
      7  903  743 (S21E01 L= 73)   237     85        0 ( 859/ 721- 936/ 769)  53.1 E03-W00 S23-S18
      8 1632  593 (S14W52 L=126)   425     46        1 (1609/ 570-1659/ 615)   0.0 W49-W55 S15-S12
      9 1243  540 (S26W26 L=100)   870    238        3 (1157/ 509-1350/ 606)  45.0 W20-W30 S29-S19
     10  971  305 (S50W18 L= 93)   313     34        1 ( 945/ 285- 993/ 327)   0.0 W15-W21 S51-S47
     11 1220  239 (S48W45 L=119)  1677    192        4 (1147/ 206-1293/ 292)   3.0 W36-W56 S50-S43
     12 1324  233 (S45W56 L=130)   540     53        1 (1307/ 210-1345/ 268) -32.0 W49-W60 S47-S42
  

    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