Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1564 1635 (N39W51 L=113)   431     85        2 (1540/1610-1586/1653) -33.7 W49-W53 N37-N40
      2 1293 1402 (N26W21 L= 83)   323     69        1 (1256/1382-1334/1415)  -9.5 W18-W24 N25-N27
      3  467 1094 (N12E37 L= 25)  1106    227        6 ( 418/1017- 524/1211) -54.0 E40-E32 N06-N18
      4 1473  718 (S18W28 L= 89)   564    167        0 (1395/ 686-1544/ 773)  21.8 W22-W32 S19-S14
      5 1591  570 (S29W39 L=100)  1209    231        5 (1516/ 499-1674/ 657)  39.3 W33-W44 S33-S23
      6 1460  435 (S38W30 L= 92)   543    117        2 (1437/ 398-1494/ 481) -66.8 W27-W33 S41-S34
      7 1394  387 (S42W25 L= 87)   691    112        1 (1343/ 364-1442/ 407) -16.7 W20-W30 S43-S39
  

    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