Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  338 1501 (N30E61 L=180)   643     99        0 ( 309/1475- 382/1533) -12.1 E68-E54 N28-N32
      2  383 1253 (N15E46 L=196)   390    130        0 ( 368/1185- 410/1296) -74.1 E47-E43 N10-N17
      3  315  650 (S24E55 L=187)   752    157        0 ( 292/ 575- 355/ 707)  42.3 E60-E48 S29-S19
      4  622  507 (S32E29 L=213)   587    192        0 ( 583/ 435- 686/ 548)  39.3 E34-E23 S38-S29
      5 1456  395 (S39W39 L=281)   503    122        0 (1397/ 371-1510/ 422)  24.0 W34-W43 S40-S36
      6  969  369 (S42E02 L=239)   613     85        2 ( 927/ 338-1006/ 395)   0.0 E05-W00 S44-S39
      7  687  358 (S44E28 L=214)   406     95        0 ( 654/ 336- 728/ 388) -38.7 E30-E24 S45-S41
      8  784  327 (S46E19 L=223)   311     46        1 ( 761/ 303- 803/ 350)  53.1 E21-E16 S48-S44
  

    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