Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1625 1363 (N17W42 L=276)  1185    326        0 (1527/1240-1726/1483) -53.1 W36-W48 N09-N25
      2  428 1413 (N23E41 L=193)   268     93        0 ( 415/1368- 446/1455)  78.7 E42-E40 N20-N25
      3 1235 1393 (N19W14 L=248)   241     58        1 (1209/1376-1264/1418) -33.7 W12-W16 N18-N21
      4 1885 1285 (N13W69 L=303)   543     57        1 (1869/1259-1907/1316)  33.7 W65-W72 N11-N15
      5  344  799 (S14E48 L=186)  1180    275        1 ( 317/ 667- 375/ 925)  72.6 E52-E45 S22-S06
      6  816  786 (S17E14 L=221)   250    129        0 ( 772/ 768- 874/ 826) -26.6 E16-E09 S17-S14
      7  698  642 (S25E23 L=211)   204     68        0 ( 668/ 618- 735/ 661)  21.8 E25-E20 S26-S24
      8 1601  594 (S30W43 L=277)   320     65        0 (1570/ 566-1637/ 622)  31.0 W40-W45 S31-S28
      9 1693  469 (S38W60 L=294)   662     65        2 (1673/ 433-1708/ 505) -29.1 W56-W65 S40-S35
     10 1429  461 (S39W31 L=266)  1042    236        2 (1335/ 433-1547/ 496)   6.3 W23-W41 S41-S36
     11  584  259 (S53E54 L=180)   539     69        2 ( 567/ 226- 610/ 289)  24.8 E64-E48 S56-S50
  

    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