Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1495 1599 (N23W50 L=180)   454     97        1 (1457/1569-1541/1623) -68.2 W48-W52 N20-N25
      2 1357 1463 (N20W33 L=163)  1205    468        2 (1161/1419-1543/1532) -34.3 W23-W45 N13-N29
      3 1795 1329 (S01W61 L=190)   562     71        1 (1762/1303-1825/1359) -36.9 W58-W63 S02-N02
      4 1121 1237 (N15W11 L=140)   290    141        0 (1053/1226-1191/1247) -26.6 W07-W15 N12-N16
      5  794 1185 (N20E09 L=120)   378    103        2 ( 740/1162- 849/1203) -23.2 E12-E05 N18-N22
      6 1896 1060 (S19W63 L=193)  1505    204        2 (1868/ 973-1917/1162)  80.5 W58-W67 S23-S11
      7  309 1089 (N26E47 L= 83)   353     64        2 ( 280/1073- 349/1110) -26.6 E49-E43 N24-N27
      8  185 1068 (N28E62 L= 67)   413     49        3 ( 162/1042- 203/1095) -53.1 E65-E59 N25-N29
      9  571  982 (N14E28 L=102)   495    294        0 ( 460/ 947- 691/1010) -17.4 E35-E19 N11-N17
     10 1289  635 (S24W05 L=134)   817    209        0 (1261/ 546-1304/ 753)  59.7 W01-W08 S29-S16
     11  510  450 (S15E51 L= 79)   841    216        1 ( 481/ 368- 552/ 538)  63.4 E54-E46 S21-S09
     12 1236  467 (S33E03 L=126)   238     36        1 (1215/ 447-1258/ 486)   0.0 E04-E01 S34-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