Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1309 1668 (N35W43 L=157)  2407    456        3 (1128/1600-1506/1704) -13.2 W26-W60 N28-N39
      2 1113 1519 (N33W21 L=135)   629    182        1 (1078/1433-1156/1598)  71.6 W18-W24 N26-N38
      3  450 1068 (N23E36 L= 79)   246     60        0 ( 423/1049- 480/1093) -36.9 E37-E33 N21-N24
      4 1927 1017 (S21W69 L=183)   698     63        1 (1904/ 988-1954/1046) -31.0 W64-W76 S24-S18
      5 1676  869 (S20W37 L=151)   263     54        0 (1658/ 841-1692/ 901)  56.3 W34-W38 S21-S18
      6 1377  852 (S12W16 L=130)   269     96        1 (1333/ 824-1424/ 879) -45.0 W13-W18 S14-S09
      7 1697  738 (S29W37 L=151)  2050    365        9 (1567/ 601-1839/ 839)  14.0 W27-W53 S36-S20
      8 1236  611 (S22W01 L=116)   652    250        0 (1164/ 513-1282/ 718)  42.3 E05-W06 S26-S16
      9  488  411 (S17E55 L= 59)   392     79        0 ( 451/ 385- 513/ 435) -53.1 E57-E53 S19-S14
     10  633  323 (S26E50 L= 65)   350     55        0 ( 611/ 300- 656/ 348)  20.6 E53-E45 S27-S24
  

    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