Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  334 1503 (N55E52 L= 13)   498     50        0 ( 315/1474- 352/1528) -39.8 E56-E47 N52-N57
      2  552 1497 (N48E22 L= 43)   341     39        1 ( 527/1470- 576/1517)  45.0 E23-E18 N46-N50
      3 1892 1235 (S08W70 L=134)   751     44        0 (1867/1211-1915/1257) -29.7 W66-W73 S10-S06
      4 1810  898 (S24W50 L=114)  3246    595        2 (1689/ 668-1877/1150)  48.8 W35-W60 S34-S09
      5  217  881 (N16E61 L=  3)  4638    532        7 ( 149/ 714- 376/1110) -55.3 E70-E49 N04-N30
      6 1200  514 (S28E04 L= 61)   249     44        0 (1173/ 496-1226/ 536)  14.0 E05-E01 S28-S26
      7 1035  338 (S35E21 L= 43)  1021    196        0 ( 965/ 278-1089/ 430)  25.0 E28-E13 S38-S29
      8  449  362 (S20E64 L=360)   737     88        3 ( 419/ 339- 480/ 388)  -8.1 E69-E59 S21-S17
      9  857  177 (S43E48 L= 16)  1266    138        1 ( 824/ 117- 895/ 229)  21.8 E64-E39 S49-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