Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1698 1520 (N20W57 L= 57)   367    111        0 (1665/1503-1726/1539)   8.1 W53-W60 N19-N21
      2 1819 1436 (N14W66 L= 67)   866    152        1 (1775/1401-1862/1475) -32.0 W62-W72 N11-N16
      3  231 1270 (N25E54 L=307)   643    170        1 ( 185/1245- 280/1287) -16.7 E58-E48 N23-N27
      4  219 1067 (N13E55 L=305)  3515    721       10 ( 103/ 910- 321/1229)  37.7 E71-E45 N04-N22
      5 1602  870 (S17W35 L= 36)   254     91        1 (1584/ 822-1618/ 914)  68.2 W33-W36 S19-S14
      6  778  722 (S15E19 L=341)   448    169        0 ( 746/ 648- 812/ 809)  52.1 E22-E15 S18-S09
      7 1634  780 (S23W37 L= 38)   252     65        0 (1609/ 762-1661/ 806) -36.9 W35-W39 S24-S20
      8  643  738 (S12E27 L=333)   290     99        0 ( 629/ 685- 655/ 787)  71.6 E29-E25 S14-S08
      9 1276  604 (S29W10 L= 11)   880    176        2 (1242/ 530-1313/ 706)  85.2 W07-W13 S33-S21
     10  422  655 (S14E45 L=315)   678    111        2 ( 382/ 595- 448/ 703) -81.9 E47-E43 S17-S10
     11  288  534 (S19E63 L=298)  1913    261        5 ( 258/ 456- 322/ 618)  51.3 E67-E57 S24-S14
     12 1053  465 (S35E07 L=354)   267     88        0 (1038/ 426-1075/ 512)  81.9 E07-E05 S38-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