Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1621 1361 (N33W30 L=124)   442     69        0 (1606/1322-1636/1400) -68.2 W28-W31 N30-N35
      2 1804 1104 (N23W51 L=145)  2531    399        5 (1688/ 976-1905/1190)  43.5 W42-W66 N12-N31
      3  395  772 (S35E37 L= 57)  2687    358        7 ( 310/ 669- 482/ 921)  52.7 E47-E31 S44-S23
      4 1754  658 (S04W61 L=154)   551     63        1 (1735/ 627-1770/ 700) -63.4 W58-W63 S05-S01
      5 1379  604 (S18W34 L=128)   638    186        1 (1308/ 569-1448/ 671)  58.0 W30-W35 S21-S12
      6 1151  606 (S24W20 L=114)   370     70        0 (1121/ 586-1192/ 628)   0.0 W17-W22 S25-S22
      7 1449  559 (S18W41 L=135)   618    121        0 (1390/ 523-1503/ 589)  50.2 W38-W43 S21-S15
      8 1018  527 (S33W15 L=109)   711    218        2 ( 944/ 471-1112/ 584)  61.4 W11-W18 S38-S26
      9  682  501 (S46E08 L= 86)   251     61        1 ( 663/ 470- 702/ 537) -80.5 E08-E06 S48-S42
     10 1636  510 (S15W59 L=152)   464     78        0 (1599/ 487-1670/ 537)  51.3 W56-W61 S17-S12
  

    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