Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  748 1503 (N42E09 L=351)  1435    318        0 ( 662/1443- 895/1595) -46.8 E16-W00 N34-N50
      2 1172 1348 (N23W16 L= 16)   203    122        0 (1158/1294-1195/1394)  71.6 W15-W18 N19-N25
      3  329 1213 (N30E47 L=312)  2602    443        8 ( 168/1086- 471/1335) -23.2 E76-E34 N20-N39
      4 1751  945 (S14W47 L= 46)   522     96        0 (1726/ 891-1776/ 989)  38.7 W43-W49 S17-S11
      5  180  936 (N14E66 L=294)   607    116        1 ( 160/ 881- 195/ 976) -68.2 E68-E63 N11-N16
      6  566  958 (N11E30 L=330)   249     75        0 ( 526/ 938- 604/ 976)   0.0 E32-E26 N10-N12
      7  250  780 (N04E60 L=300)   434     64        3 ( 220/ 758- 273/ 803)  18.4 E64-E57 N03-N06
      8  794  676 (S09E21 L=339)   251     67        0 ( 765/ 656- 824/ 694) -45.0 E22-E19 S10-S07
      9  721  613 (S12E27 L=333)   284    138        1 ( 668/ 573- 776/ 658)  18.4 E31-E22 S13-S09
     10 1264  429 (S36W03 L=  3)   482    105        1 (1217/ 408-1307/ 452)   0.0 W00-W06 S37-S34
     11  934  348 (S34E23 L=337)   597    139        2 ( 911/ 285- 959/ 416)  48.0 E27-E18 S39-S29
  

    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