Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  551 1673 (N20E51 L=268)  1018    156        1 ( 466/1635- 602/1721)  53.1 E57-E48 N16-N24
      2  811 1668 (N25E32 L=287)   622     86        2 ( 782/1634- 853/1712) -23.2 E35-E28 N23-N28
      3 1225 1362 (N17W03 L=322)   367     82        1 (1200/1321-1252/1397)  90.0 W00-W04 N13-N19
      4 1833 1172 (N25W52 L= 11)  1572    299        1 (1783/1028-1873/1293) -47.3 W45-W58 N17-N31
      5 1113  872 (S13W09 L=329)   225     57        0 (1092/ 847-1146/ 894)  56.3 W07-W10 S14-S11
      6  907  641 (S32W04 L=323)   231     83        0 ( 868/ 617- 942/ 662)  51.3 W01-W05 S34-S29
      7 1133  533 (S32W22 L=342)   510    137        0 (1091/ 502-1185/ 590) -11.3 W16-W26 S33-S29
      8  605  553 (S46E16 L=303)   282     42        0 ( 579/ 539- 643/ 572)   0.0 E19-E12 S47-S45
      9  636  352 (S61E10 L=310)   355     37        0 ( 616/ 328- 655/ 376)  68.2 E12-E07 S63-S58
     10  796  283 (S61W12 L=332)   502     72        1 ( 765/ 255- 834/ 315)  -5.7 W08-W18 S62-S57
  

    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