Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  430 1703 (N32E57 L=202)   617    167        0 ( 413/1682- 457/1737) -29.1 E61-E52 N30-N35
      2 1119 1471 (N22W00 L=259)   391     87        1 (1085/1435-1165/1496)  36.9 E01-W03 N20-N24
      3  203 1360 (N08E66 L=193)  2065    331        4 ( 131/1317- 261/1393)   9.5 E76-E58 N05-N09
      4 1386 1227 (N11W19 L=278)  1203    220        3 (1329/1126-1433/1310) -55.0 W15-W23 N05-N15
      5  247 1178 (S03E55 L=203)   518    123        1 ( 219/1129- 284/1224) -56.3 E58-E51 S06-N00
      6  158 1088 (S09E64 L=195)   635    105        2 ( 136/1049- 189/1122)  45.0 E66-E60 S11-S06
      7 1838 1076 (N09W56 L=314)   372     82        2 (1821/1038-1856/1113)  78.7 W54-W57 N06-N11
      8  594  942 (S15E26 L=233)   261    103        0 ( 542/ 903- 620/ 981) -38.7 E29-E23 S17-S13
      9  693  842 (S20E19 L=240)   197     47        1 ( 673/ 818- 719/ 867) -45.0 E19-E16 S21-S18
     10 1487  806 (S12W32 L=290)   418     97        2 (1433/ 782-1537/ 833)  33.7 W28-W35 S14-S10
     11 1344  758 (S17W23 L=282)   217     45        0 (1312/ 739-1372/ 776)   0.0 W21-W25 S18-S16
     12  866  203 (S66W01 L=260)   652     87        1 ( 840/ 165- 892/ 227) -41.2 E03-W05 S70-S63
  

    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