Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1038 1652 (N39E14 L=194)   608     97        1 (1009/1610-1075/1697) -32.0 E17-E09 N36-N42
      2 1341 1475 (N33W12 L=220)   955    168        2 (1276/1412-1405/1529) -33.7 W07-W17 N29-N36
      3 1558 1186 (N19W31 L=239)   474    158        0 (1532/1114-1598/1258) -53.1 W28-W35 N14-N22
      4 1332 1116 (N10W17 L=225)   380     88        2 (1301/1088-1356/1154)  63.4 W15-W18 N08-N12
      5 1837 1026 (N15W59 L=267)   782     78        3 (1805/ 991-1864/1056) -23.2 W55-W63 N12-N16
      6  369  893 (S21E42 L=166)  2353    443        3 ( 217/ 755- 503/1028)  38.2 E61-E32 S32-S10
      7 1261  675 (S17W22 L=230)   338     38        0 (1235/ 651-1287/ 698) -14.0 W19-W23 S18-S16
      8 1733  625 (S11W60 L=268)  1556    159        3 (1650/ 591-1791/ 666)   6.3 W52-W70 S14-S08
      9 1522  495 (S23W47 L=255)   349     53        1 (1497/ 471-1546/ 519)  56.3 W45-W48 S24-S21
     10 1162  266 (S48W35 L=243)   372     67        0 (1139/ 245-1202/ 292) -12.5 W31-W40 S49-S46
     11  950  261 (S55W17 L=225)   699     80        1 ( 917/ 240- 989/ 287)  -6.3 W12-W21 S56-S52
     12  749  262 (S61E05 L=203)   510     85        1 ( 709/ 249- 790/ 284)   5.7 E10-W00 S63-S58
  

    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