Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1236 1392 (N20W23 L= 98)   467    132        0 (1216/1331-1261/1455)  40.6 W20-W27 N17-N23
      2 1803 1261 (S04W58 L=133)  1865    239        3 (1763/1162-1845/1382)  76.0 W54-W61 S10-N02
      3 1029 1185 (N15W05 L= 80)   442    174        1 ( 945/1165-1123/1216) -16.7 W00-W10 N12-N17
      4  633 1051 (N18E22 L= 53)   696    229        0 ( 541/1025- 751/1086) -30.3 E27-E15 N13-N21
      5 1675 1018 (S13W39 L=114)   302     79        1 (1657/ 976-1691/1056)  53.1 W37-W40 S15-S11
      6 1578  848 (S19W28 L=103)  1432    338        4 (1536/ 720-1604/1043)  59.5 W22-W34 S26-S09
      7 1141  654 (S18E04 L= 71)   311    104        2 (1095/ 628-1179/ 674)   0.0 E06-E01 S18-S16
      8 1319  560 (S29W04 L= 79)  1222    272        2 (1204/ 494-1387/ 670)  17.1 E03-W11 S32-S23
      9  759  594 (S11E27 L= 48)   239     94        0 ( 715/ 575- 808/ 623)   0.0 E30-E23 S11-S10
     10  861  519 (S18E24 L= 51)   378    151        0 ( 837/ 459- 885/ 590)  45.0 E27-E19 S21-S14
     11  532  416 (S17E51 L= 24)   372     75        2 ( 496/ 397- 568/ 435) -45.0 E53-E48 S18-S14
     12  564  327 (S23E54 L= 21)   584    140        2 ( 523/ 300- 610/ 354) -15.9 E57-E50 S24-S21
  

    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