Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1404 1561 (N32W38 L= 32)   464    171        2 (1327/1539-1470/1589) -10.3 W32-W44 N30-N33
      2 1099 1435 (N30W12 L=  7)   812    303        2 (1035/1382-1186/1477)   9.5 W07-W19 N27-N31
      3 1667 1354 (N14W50 L= 45)  1280    214        3 (1595/1289-1739/1398) -48.0 W44-W54 N08-N18
      4 1135  962 (N01W06 L=  0)  1754    467        7 (1035/ 779-1185/1173) -87.6 W02-W07 S10-N14
      5  523 1070 (N16E32 L=323)   552    142        1 ( 462/1036- 581/1112) -26.6 E35-E27 N13-N18
      6 1441  901 (S08W24 L= 18)   786    262        3 (1370/ 829-1542/ 956) -47.7 W19-W29 S14-S03
      7 1310  795 (S12W14 L=  8)   385    104        0 (1290/ 741-1333/ 849)  74.1 W12-W15 S15-S08
      8 1317  582 (S26W12 L=  6)   634    116        1 (1285/ 522-1348/ 639)  54.5 W09-W14 S29-S22
      9  968  347 (S37E17 L=337)  1187    297        1 ( 889/ 280-1016/ 460)  45.0 E25-E11 S42-S28
  

    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