Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  365 1611 (N29E56 L=322)   869    153        1 ( 345/1557- 380/1667) -71.6 E59-E53 N24-N33
      2 1514 1533 (N32W30 L= 48)   528     81        1 (1487/1507-1539/1566)  45.0 W28-W32 N29-N34
      3  355 1490 (N20E52 L=326)   438     98        1 ( 322/1470- 390/1511)  18.4 E54-E48 N18-N21
      4 1256 1356 (N17W12 L= 30)   505    133        2 (1224/1290-1284/1431) -77.5 W09-W13 N12-N21
      5  646 1339 (N11E26 L=352)   880    243        1 ( 572/1244- 765/1397) -29.7 E31-E17 N06-N14
      6  246 1247 (N05E56 L=322)  1222    210        4 ( 212/1175- 295/1298)  40.6 E60-E51 N00-N07
      7 1509 1182 (N09W29 L= 47)   558    128        0 (1449/1119-1559/1228) -32.0 W24-W32 N05-N11
      8 1919 1078 (N08W67 L= 85)   692    129        3 (1898/1043-1939/1106) -45.0 W63-W70 N06-N09
      9 1786  982 (N00W52 L= 70)   689    147        0 (1764/ 909-1804/1058) -54.5 W49-W54 S03-N04
     10  301  749 (S25E50 L=328)  2010    479        1 ( 228/ 594- 402/ 916)  48.0 E60-E40 S35-S15
     11  480  443 (S46E42 L=336)   352     76        0 ( 461/ 415- 498/ 475)  38.7 E44-E38 S47-S43
  

    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