Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  963 1872 (N45E38 L=  8)  1095    110        2 ( 921/1843-1016/1899)  38.7 E41-E33 N42-N49
      2 1063 1696 (N34E19 L= 27)  2306    512        6 (1003/1543-1103/1857) -45.0 E32-E10 N24-N46
      3  383 1655 (N15E66 L=340)   724    120        3 ( 356/1624- 411/1689)  78.7 E69-E63 N12-N17
      4  345 1492 (N04E59 L=347)   529    156        1 ( 291/1457- 398/1520)  32.0 E63-E54 N01-N06
      5  674 1413 (N06E31 L= 15)   804    201        0 ( 610/1345- 739/1486) -22.6 E37-E24 N03-N09
      6 1026 1322 (N09E08 L= 38)   371    113        0 ( 963/1287-1079/1346)   0.0 E11-E04 N08-N10
      7  496 1317 (S03E39 L=  7)   310     83        0 ( 460/1295- 531/1335)   0.0 E42-E36 S04-S02
      8  978 1091 (S04E04 L= 42)   905    284        0 ( 847/ 969-1083/1193)  67.6 E08-W00 S14-N03
      9  130 1020 (S27E65 L=341)   736    164        3 ( 107/ 970- 159/1069)  90.0 E67-E62 S30-S23
     10  544 1009 (S20E28 L= 18)   237     57        0 ( 512/ 989- 575/1028)   0.0 E30-E25 S20-S18
     11 1876  921 (N14W60 L=106)   748    203        4 (1849/ 896-1913/ 967)  -6.3 W56-W65 N11-N16
     12  519  918 (S26E28 L= 18)   327     35        1 ( 498/ 891- 542/ 945)  76.0 E29-E26 S27-S23
  

    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