Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  705 1305 (N24E20 L=195)   242     68        1 ( 679/1272- 727/1335) -53.1 E22-E18 N21-N25
      2  771 1171 (N15E15 L=200)   460    159        3 ( 699/1142- 850/1199) -11.3 E20-E10 N13-N17
      3  693  937 (N01E21 L=194)   244     48        1 ( 665/ 922- 726/ 956)   0.0 E23-E19 N00-N02
      4 1191  895 (S04W09 L=225)   671    217        2 (1090/ 861-1285/ 934) -22.6 W03-W15 S06-S01
      5  194  777 (S09E67 L=148)   918     98        1 ( 169/ 744- 218/ 821) -63.4 E70-E62 S10-S06
      6 1621  670 (S22W41 L=256)   498    130        0 (1580/ 621-1654/ 732)  66.8 W37-W43 S25-S18
      7  334  636 (S18E54 L=161)   581    122        0 ( 307/ 587- 360/ 691) -81.9 E56-E52 S20-S13
      8  591  569 (S22E33 L=182)  1310    273        3 ( 523/ 469- 687/ 654)  33.7 E42-E24 S29-S16
      9 1615  490 (S36W45 L=261)   500     42        1 (1584/ 465-1644/ 508)   0.0 W42-W48 S37-S34
     10  539  380 (S37E46 L=169)   797    120        1 ( 522/ 318- 570/ 434)  51.3 E53-E40 S42-S32
  

    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