Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1549 1697 (N26W63 L= 48)   781     91        1 (1494/1676-1593/1712) -24.0 W58-W68 N24-N28
      2 1315 1593 (N28W37 L= 22)   269     61        1 (1290/1562-1341/1618) -90.0 W35-W38 N25-N30
      3 1631 1507 (N14W55 L= 40)  1135    288        2 (1506/1452-1755/1536) -33.7 W46-W64 N07-N19
      4 1139 1408 (N23W18 L=  3)   726    262        0 (1016/1362-1245/1441) -42.3 W11-W22 N17-N28
      5 1629 1293 (N03W44 L= 29)   295     68        1 (1599/1260-1656/1319)  18.4 W41-W47 N01-N04
      6  175 1006 (N23E63 L=282)  1225    156        3 ( 110/ 974- 255/1031)  -2.5 E75-E52 N21-N25
      7 1394 1002 (S06W20 L=  5)   267     46        1 (1365/ 982-1416/1030)  33.7 W18-W21 S07-S04
      8 1308  847 (S12W11 L=356)   223     72        0 (1281/ 826-1331/ 875)  18.4 W09-W12 S13-S10
      9 1491  626 (S31W18 L=  3)  1766    286        8 (1443/ 512-1540/ 738)  50.2 W12-W23 S38-S23
     10 1598  550 (S41W26 L= 11)   475     71        0 (1563/ 524-1631/ 572)   9.5 W23-W29 S42-S39
     11  931  327 (S34E29 L=316)   395     78        0 ( 917/ 290- 951/ 368)  56.3 E30-E25 S37-S31
  

    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