Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1075 1341 (N21W12 L= 78)   214    103        0 (1039/1311-1122/1358)   0.0 W08-W14 N20-N22
      2  477 1227 (N31E30 L= 36)  1886    357        2 ( 345/1179- 633/1268) -14.0 E42-E18 N24-N34
      3  404 1011 (N20E39 L= 27)  1202    307        1 ( 285/ 957- 531/1086) -30.1 E49-E30 N14-N26
      4  679  930 (N09E22 L= 44)  1097    328        0 ( 630/ 868- 752/1031)  71.6 E24-E16 N04-N14
      5  177  943 (N20E64 L=  2)  1373    189        4 ( 138/ 860- 237/1023) -47.7 E69-E57 N14-N25
      6 1814  892 (S25W50 L=116)   385     60        0 (1788/ 870-1836/ 920)   9.5 W46-W52 S26-S24
      7 1677  561 (S43W34 L=100)  1230    231        3 (1576/ 505-1776/ 619)  -2.5 W23-W46 S44-S40
      8  435  380 (S18E63 L=  3)  3284    383        4 ( 315/ 250- 539/ 497)  78.7 E71-E57 S28-S08
      9  912  276 (S37E34 L= 33)  1842    230        4 ( 824/ 199- 986/ 354)  13.4 E45-E23 S41-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