Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1028 1869 (N58W46 L=328)  1836    257        3 ( 953/1782-1113/1955)   6.6 W29-W81 N52-N60
      2  629 1636 (N55E10 L=272)   609    106        2 ( 581/1601- 670/1669)  11.3 E15-E05 N53-N57
      3 1598 1509 (N15W53 L=335)   393     67        1 (1572/1487-1633/1535) -26.6 W50-W55 N13-N17
      4  584 1411 (N40E19 L=264)   330     55        1 ( 549/1387- 616/1429)   0.0 E21-E15 N38-N41
      5 1778  960 (S20W47 L=330)   582    184        0 (1741/ 882-1818/1036)  35.0 W42-W52 S23-S16
      6 1845  959 (S22W54 L=337)   479    103        0 (1816/ 916-1872/1003)  20.6 W50-W58 S24-S20
      7 1468  708 (S25W18 L=301)   768    289        1 (1336/ 642-1581/ 757)  -6.0 W07-W26 S27-S22
      8 1706  648 (S37W37 L=320)   341     38        1 (1682/ 625-1725/ 672) -76.0 W34-W39 S39-S35
      9 1247  466 (S33E03 L=280)   583    148        1 (1170/ 448-1306/ 497) -18.4 E07-W02 S35-S30
     10 1601  436 (S50W26 L=309)   901     88        0 (1566/ 390-1631/ 475) -55.0 W22-W30 S55-S45
     11  879  398 (S27E28 L=254)   565    117        1 ( 861/ 346- 907/ 454)  36.9 E31-E23 S30-S23
  

    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