Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1386 1706 (N43W48 L= 13)   408     89        1 (1356/1682-1417/1725)   0.0 W42-W51 N41-N44
      2  612 1432 (N37E25 L=300)   285     85        1 ( 579/1399- 649/1461) -45.0 E27-E21 N34-N39
      3 1539 1287 (N14W38 L=  3)   822    181        1 (1508/1207-1584/1362) -90.0 W35-W39 N08-N18
      4 1148 1219 (N16W10 L=336)   491    192        2 (1058/1172-1228/1259) -37.9 W05-W14 N12-N19
      5  402 1162 (N22E41 L=284)   470    136        2 ( 360/1117- 456/1207) -33.7 E45-E36 N18-N24
      6  812 1147 (N16E12 L=314)   283     56        0 ( 792/1115- 832/1176)  76.0 E12-E10 N14-N18
      7 1734  990 (S08W48 L= 13)   656    156        3 (1690/ 940-1779/1048)  24.0 W43-W52 S10-S04
      8 1769  941 (S11W51 L= 16)   367     87        1 (1739/ 906-1788/ 979)  31.0 W47-W52 S13-S09
      9 1687  709 (S25W43 L=  8)   887    206        2 (1630/ 624-1742/ 792)  45.0 W37-W47 S29-S19
     10 1554  749 (S19W31 L=356)   275     70        0 (1542/ 707-1568/ 790)  78.7 W30-W32 S21-S16
     11  486  401 (S28E54 L=271)  2428    230        6 ( 455/ 288- 533/ 500)  29.4 E76-E44 S38-S20
     12  957  323 (S40E18 L=307)   300     55        0 ( 928/ 306- 989/ 342)  18.4 E20-E14 S41-S38
  

    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