Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  798 1398 (N34E03 L= 96)   330     72        1 ( 761/1380- 842/1415) -11.3 E05-W00 N32-N34
      2  895 1207 (N20E02 L= 96)   309    108        2 ( 873/1147- 918/1258)  56.3 E04-W00 N16-N22
      3 1881 1131 (S14W61 L=160)  2077    359        4 (1838/1024-1908/1247)  47.7 W53-W68 S19-S08
      4 1549 1119 (S03W33 L=132)   935    214        2 (1521/1019-1590/1227)  74.7 W30-W36 S09-N02
      5  323  945 (N19E47 L= 52)   388    173        0 ( 260/ 909- 376/ 992) -41.2 E51-E43 N15-N22
      6 1768  861 (S25W43 L=142)   720    136        0 (1731/ 800-1805/ 936)  32.0 W39-W47 S28-S21
      7 1239  659 (S20W02 L=101)  1055    392        2 (1177/ 522-1279/ 881)  52.1 E05-W10 S27-S08
      8  731  818 (N02E22 L= 77)   212     57        1 ( 714/ 791- 759/ 840)  18.4 E23-E19 N00-N03
      9 1034  645 (S15E10 L= 89)   391    187        0 (1019/ 556-1070/ 724)  52.1 E12-E05 S20-S11
     10  346  618 (S01E56 L= 43)   550    102        2 ( 332/ 573- 364/ 665)  38.7 E58-E52 S03-N02
     11  990  397 (S29E21 L= 77)  1022    228        1 ( 897/ 356-1058/ 474)  -4.8 E26-E14 S32-S25
  

    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