Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1383 1707 (N33W47 L=137)   320     80        0 (1344/1691-1413/1726)   0.0 W43-W49 N32-N34
      2 1871 1229 (S07W64 L=155)  1331    261        3 (1825/1136-1914/1353) -85.6 W61-W67 S13-N00
      3 1208 1281 (N14W17 L=107)   225     80        0 (1165/1246-1252/1308) -63.4 W15-W18 N10-N16
      4 1714 1167 (S06W46 L=137)  1155    259        2 (1683/1035-1735/1302)  70.3 W43-W49 S13-N01
      5  997 1158 (N12W02 L= 92)   295     78        0 ( 947/1142-1044/1178) -21.8 E01-W04 N11-N13
      6  401 1060 (N20E38 L= 52)   463    121        0 ( 349/1044- 457/1078) -12.5 E42-E33 N18-N21
      7  571 1023 (N14E27 L= 64)   303     88        0 ( 533/ 979- 597/1057) -63.4 E28-E25 N11-N17
      8  137  911 (N17E67 L= 23)  1540    241        3 ( 101/ 857- 170/ 969) -60.3 E72-E62 N13-N20
      9 1126  647 (S20E04 L= 87)   580    205        2 (1095/ 553-1149/ 746)  55.0 E07-W00 S25-S14
     10 1564  671 (S31W25 L=115)   770    119        3 (1537/ 619-1593/ 726)  81.9 W22-W27 S33-S26
     11  510  637 (S06E41 L= 49)   848    177        1 ( 460/ 597- 598/ 704) -53.1 E44-E36 S09-S01
     12 1309  572 (S30W06 L= 96)   400    107        2 (1287/ 522-1329/ 636)  60.3 W03-W08 S32-S25
  

    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