Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1506 1607 (N27W45 L=239)  1131    207        2 (1433/1571-1565/1640)   4.4 W38-W51 N25-N28
      2  852 1563 (N37E01 L=192)   816    128        2 ( 790/1508- 915/1591) -14.0 E05-W03 N33-N39
      3 1356 1561 (N27W33 L=226)   271     52        1 (1343/1528-1374/1584)  36.9 W30-W34 N25-N28
      4 1609 1358 (N10W44 L=237)   352     78        0 (1577/1333-1651/1384) -38.7 W41-W46 N08-N12
      5  599 1277 (N23E23 L=171)   788    183        1 ( 517/1227- 671/1325) -42.0 E28-E17 N18-N27
      6 1296 1270 (N11W20 L=213)   249     69        1 (1260/1244-1331/1299) -53.1 W18-W22 N08-N12
      7  924 1235 (N15E02 L=191)   887    258        0 ( 852/1142- 989/1289)  41.2 E08-W01 N10-N17
      8  218 1224 (N26E54 L=139)   533    100        0 ( 192/1199- 261/1245)   8.1 E56-E49 N24-N27
      9  782  424 (S31E29 L=165)   463    105        0 ( 735/ 396- 823/ 456)  14.0 E32-E24 S32-S29
     10  961  210 (S52E28 L=165)   368     36        0 ( 935/ 193- 985/ 227) -45.0 E30-E26 S54-S50
  

    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