Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1339 1574 (N33W12 L=356)   387     98        1 (1285/1560-1387/1591)  15.9 W08-W15 N31-N34
      2 1180 1407 (N19W04 L=348)   541    149        1 (1112/1380-1259/1429)   7.1 W00-W08 N17-N19
      3 1260 1303 (N13W10 L=354)  2374    277        8 (1108/1202-1385/1382)   0.0 W00-W17 N06-N19
      4 1009 1297 (N09E05 L=339)   374     76        1 ( 975/1261-1038/1323) -45.0 E06-E03 N07-N10
      5 1538  957 (S02W32 L= 16)   378    113        0 (1474/ 942-1591/ 971)  23.2 W27-W35 S03-N00
      6 1534  548 (S25W44 L= 28)  4514   1035        6 (1224/ 403-1837/ 965)  52.2 W23-W58 S37-N04
      7 1224  901 (S10W13 L=357)   196     65        0 (1205/ 872-1238/ 942)  90.0 W12-W14 S12-S07
      8  315  853 (S25E47 L=297)  1340    289        3 ( 230/ 778- 449/ 904)  20.0 E57-E35 S29-S20
      9 1057  706 (S24W07 L=351)  3068    995        4 ( 752/ 632-1508/ 858)   3.5 E13-W36 S29-S17
     10  416  547 (S44E43 L=301)  1623    356        3 ( 262/ 516- 596/ 568)   2.8 E65-E24 S45-S40
     11 1770  527 (S21W69 L= 52)   544     89        0 (1744/ 486-1791/ 567) -78.7 W66-W70 S23-S18
  

    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