Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  269 1480 (N52E59 L=188)   688     77        0 ( 232/1463- 313/1491)  -3.2 E69-E50 N50-N53
      2  698 1344 (N31E12 L=235)   238     64        0 ( 661/1331- 733/1357) -11.3 E14-E09 N29-N31
      3  424 1243 (N31E35 L=213)   436     89        0 ( 373/1212- 465/1269)   0.0 E39-E30 N30-N32
      4  455  933 (N13E37 L=210)   307     98        0 ( 414/ 912- 486/ 963)  -9.5 E39-E33 N11-N14
      5 1575  852 (S20W29 L=276)   318     95        1 (1544/ 820-1605/ 873)  11.3 W26-W31 S21-S18
      6 1716  768 (S29W39 L=287)   420     41        1 (1693/ 742-1739/ 795) -51.3 W37-W41 S31-S26
      7 1424  539 (S35W12 L=260)   282     44        0 (1406/ 515-1440/ 570)  76.0 W10-W13 S36-S32
      8 1507  313 (S58W18 L=265)   585     89        0 (1458/ 291-1546/ 336)  -4.8 W11-W23 S60-S56
      9 1274  210 (S58E12 L=235)  2299    273        3 (1139/ 175-1395/ 250) -19.8 E24-W02 S62-S52
  

    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