Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1360 1628 (N42W33 L=290)   343     83        0 (1320/1610-1397/1642) -18.4 W29-W36 N40-N43
      2 1770 1252 (N13W57 L=314)  1167    180        3 (1708/1173-1824/1306) -71.6 W51-W62 N08-N17
      3  302 1164 (N14E52 L=205)   347     53        0 ( 270/1148- 332/1179)   0.0 E55-E48 N13-N15
      4 1038  872 (S05W00 L=257)   345    147        0 ( 957/ 853-1096/ 887) -15.9 E04-W03 S06-S03
      5 1051  635 (S20W00 L=257)   288    170        1 (1018/ 572-1085/ 681) -90.0 E01-W02 S24-S17
      6 1413  639 (S22W25 L=281)   459     80        2 (1394/ 598-1443/ 673)  90.0 W23-W26 S24-S19
      7  557  593 (S22E35 L=222)   400     81        1 ( 525/ 561- 591/ 627) -45.0 E36-E32 S23-S19
      8 1306  544 (S28W18 L=274)   530     73        1 (1271/ 508-1345/ 579)  31.0 W15-W20 S30-S25
      9  535  518 (S27E39 L=218)   308    105        0 ( 513/ 473- 571/ 554)  40.6 E41-E34 S30-S24
     10  924  499 (S29E10 L=247)   693    257        0 ( 828/ 445-1038/ 551)  23.6 E17-E01 S33-S26
     11 1560  393 (S43W46 L=303)   607     47        1 (1533/ 371-1586/ 419) -24.0 W41-W50 S44-S40
  

    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